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Copyright © 1975, Texas Instruments Incorporated 


The program material contained herein is supplied without 
representation or warranty of any kind. Texas Instruments therefore 
assumes no responsibility and shall have no liability, consequential or 
otherwise, of any kind arising from use of this program material or any 
part thereof. 


ADDENDUM 


Dear Customer: 
Please mark the following corrections in the Statistics Library Program 
Manual ST1: 


PROGRAM PAGE SHOWN AS CHANGE TO 
ST 1-03-1 13 196 47 *1’ 196 87") 
ST 1-05 20 218. 18 C’ 218 18 *C’ 
ST 1-06 29 120 04 4 129 04 4 
ST 1-15 71 132 04 4 133 04 4 
ST 1-25 103 002 02 INV 002 22 INV 
ST 1-14 66-67 (see other side) 


DO NOT RERECORD THE MAGNETIC CARDS. The prerecorded magnetic 
cards have been correctly recorded at the factory. 
Change the first equation on page 18 to read: 
U;,, = 79(U; + 10°°)(Mod 10'°) 
On page 65, the minus sign should be deleted from the display column 
in Step 6. 
Note 2 on page 68 should be changed to read |X| < 21.097 
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ADDENDUM 
SR-52 Program Manual ST1 Statistics Library 


This manual was printed with an incorrect program listing for program 
ST1-14 Histogram Construction. DO NOT rerecord the magnetic card for 
this program. The correct program listing follows. 


Shares PROGRAM LISTING 
Toray [wer | oom [ ver | oom | ver [ oun [| eer 


56 
46 
14 
50 
01 
46 
19 
42 
01 
06 
94 
85 
43 
01 
04 
95 
23 


65 
43 
01 
06 
75 
43 
01 
05 
95 
80 *if pos 





“Denotes 2nd function key 


PROGRAM LISTING 


ST1-14 


“Denotes 2nd function key 





140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
tot 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 


86 *rset 
46 *LBL 
10 
43 
01 
03 
65 
43 
00 
00 
85 
43 
01 
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Table of Contents 
ee ee Toe CERO oe a ooh egal bin pat . : 
ENTERINGAPROGRAM............... iideneaiie 


ST1-01. BASIC STATISTICS FOR ONE OR TWO VARIABLES. ... . 
Calculates the means, standard deviations, and standard 
errors for one or two variables using both the normalized 

_ and unnormalized methods. Also calculates the covariance 
~ and correlation coefficient for the two variables. 


S$T1-02. PERMUTATIONS AND COMBINATIONS .......... 
Calculates the permutation and combination for 
agivennandr. 


S$T1-03-1,2. MEANS ANDMOMENTS .........+ +442. 
The arithmetic, geometric, harmonic, and generalized 
means, the first four moments, and the kurtosis and 
skewness of distribution are calculated for grouped or 
ungrouped data. 


ST1-04. RANDOM NUMBERGENERATOR............. 
Uniformly distributed and normally distributed 
random numbers are generated. | 


ST1-05. ONE WAY ANALYSIS OF VARIANCE. ........... 
Performs a one way analysis of variance on k treatment 
groups. Includes the F statistic, degrees of freedom, 
and sums of squares. 


S$T1-06. t STATISTIC EVALUATION (PAIRED OBSERVATIONS) 
. Uses t statistics to test the difference between the means 
of two normally distributed populations. 


ST1-07.t STATISTIC EVALUATION (TWO SAMPLETEST)...... 
Evaluates the t statistic with n, + n, — 2 degrees 
of freedom. 


ST el EsIGAR REGRESOION 66.6 alk whee eae ee eo 
The slope, intercept, estimated values for x or y, and 
coefficient of determination for the linear least-squares 
fit of given points (x,y) are calculated. 
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ST1-09. POWERCURVEFIT. ............80 80884 5 


The slope, intercept, estimated values for x or y, and 
correlation coefficient for a power least-squares 
fit of given points (x,y) are calculated. 


311-10. EXPONENTIAL CURVEFIT .. 22) Se a ae 


The slope, intercept, estimated values for x ory, and 
correlation coefficient for an exponential least-squares 
fit of given points (x,y) are calculated. 


ST1-11. LOGARITHMIC CURVE FIT .......... ae es 


The slope, intercept, estimated values for x or y, and 
correlation coefficient for a logarithmic least-squares 
fit of given points (x,y) are calculated. 


ST1-12-1,2. MULTIPLE LINEAR REGRESSION .......... 


The coefficients of the linear least-squares fit of 
given points (x,y,z) are calculated as well as the z 
corresponding to given x and y for this line. 


ST1-13-1,2. QUADRATIC REGRESSION ............. 


The coefficients of the quadratic least-squares fit of 
given points (x,y) are calculated and the value of y 
for a given x is determined for this curve. 


ST1-14. HISTOGRAM CONSTRUCTION ............. 


Constructs a histogram over a given range and given 
observed data points. The frequencies are calculated for 
each of 12 cells and the mean and variance are calculated 
for the entire range. 


ST1-15. NORMAL DISTRIBUTION. ........2.2.2.2.2.2.4.. | 


The standard normal distribution for a given x is 
calculated using a probability function. 


ST1-16. CHI-SQUARE DISTRIBUTION. ............. 


The chi-square density function and probability I'(4) 
are calculated, given the degrees of freedom and x. 


OU O02 C OOi ee PERPOORS EERE oe Vas eo ws ew ee ee tem ewe ws 


The integral of the t distribution is calculated, 
given the degrees of freedom and x. 
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1 4216. F DISTRIBUTION fy aaa ae ON BR aI ee eats . 80 
The probability of an event X > x (area of the “upper tail’’) 
for the F distribution is calculated for two given degrees 
of freedom and x. 


ST1-19. BIVARIATE NORMAL DISTRIBUTION ........... 84 
The bivariate density function is calculated given the 
means of x and y, standard deviations of x andy, and 
the correlation coefficient between x and y. 


$T1-20. LOGARITHMIC NORMAL DISTRIBUTION ......... 88 
The median, mode, mean, and variance are calculated 
for the logarithmic normal distribution given m and o”. 
Given an x, f(x) is then calculated for this distribution. 


oS! 21. WEIGULL DISTRIBUTION 22 4i6 cc ee eee we eke 4 90 
The probability density and cumulative probability 
functions are calculated for the Weibull distribution of 
given m,n, and x. Given Q(x), the program also calculates x. 


ST1-22. POISSON DISTRIBUTION. fa clipe Soe we we as 92 
Calculates the density function and cumulative 
distribution of the Poisson distribution given m and n. 


511-23. BINOMIAL DISTRIBUTION 2.3. 2 6s woe ee eee as 94 
The binomial density function, the cumulative distribution, 
mean, variance, and probability of at least k successes are 
calculated given n, p, and k. 


$T1-24. NEGATIVE BINOMIAL DISTRIBUTION. .......... 98 
The mean, variance, negative binomial density function, 
cumulative distribution, and probability of k or less 
successes are calculated given r, p, and k. 


S$T1-25.GEOMETRIC DISTRIBUTION .........4.4.4.2.. 102 
The variance, probability, cumulative distribution, and 
geometric density function are calculated given yu and x. 


S$T1-26. HYPERGEOMETRIC DISTRIBUTION. .......... 104 
The hypergeometric density function, cumulative 
distribution, mean, and variance are calculated 
given a, b, n, and x. 


PREFACE 


The Statistics Library includes this program Manual and 29 prerecorded 
magnetic cards. The programs in this library have been selected because 
of their common use by statisticians and other professionals who need 
rapid solutions to problems involving curve fits, distribution functions 
and general statistics. 


Each program in the Statistics Library is identified by a program name 
and number which also appears on each magnetic card. For example, 
BASIC STATISTICS FOR ONE OR TWO VARIABLES is the name of the 
first program and the number is ST1-01... where ST1 identifies the 
Statistics Library and -01 identifies the first program in the Statistics 
Library. Some programs are in two parts and have two magnetic cards 
which must be run one at atime to obtain different but related 
information. ST1-03-1 identifies the first card of the third program in the 
Statistics Library which is MEANS AND MOMENTS (1). The second card, 
MEANS AND MOMENTS (2), is identified by ST1-03-2. 


The Statistics Library programs are individually described in this manual. 
The general description of a program begins on a left-hand page, 
followed by sample problems, a pictorial of the applicable magnetic 
card, and complete user instructions. The user instructions list all the 
necessary steps to perform the specific program. The first step of the user 
instructions is enter program which should be performed according 

to the procedure on page one of this manual. Each program section also 
includes a program listing and identification of memory registers and 
flags used by the program. The program listing provides a backup 
source to check that the prerecorded program has been correctly read 
by the calculator. 


If you have difficulties in running a program: 
1. Be sure that the user instructions have been followed. 
2. Refer to Maintenance and Service Information in the Owner's 
Manual to determine the type of trouble you are having. 
3. Write the Consumer Relations Department at: 
Texas Instruments Incorporated 
P.O. Box 22283 
Dallas, Texas 75222 
or call Consumer Relations at 800-527-4980 (toll-free within all 
contiguous United States except Texas) or 800-492-4298 
(toll-free within Texas). If outside contiguous United States, call 
214-238-5461 (We regret that we cannot accept collect calls at 
this number). 


ENTERING A PROGRAM 


Use the following procedure each time you need to store the contents 
of a prerecorded magnetic card in the calculator. Use of this procedure 
is implied by the Enter program instruction in the user instructions of 
each program. 


1. Select the proper magnetic card for the program to be run. 


2. Read side A of the prerecorded 
card as follows: 
Ig (Insert card @A® ) 
Do not restrict or hold the card 
after it is caught by the drive 
motor. The display is blank until 
the calculator has completed 
reading side A. 


MMi > 
‘ Wf ., 


3. Remove the card from the left 
side of the calculator and read 
side B of the prerecorded 
card as follows: 
fy (Insert Card @#Be ) 
The disptay is blank until the 
calculator has completed 
reading side B. 


4. Remove card from the left side 
of the calculator and insert into 
upper slot such that the A side 
of the card shows in the window 
above keys A through E. 





5. |f the display flashes immediately following step 2 or 3, repeat 
the procedure beginning with step 2. 


CAUTION: Prerecorded magnetic cards may be damaged or altered if 
exposed to dust or foreign materials, permanent magnets, or 
electromagnetic fields (electric motors, power transformers, etc.). 


311-01 


BASIC STATISTICS FOR ONE OR TWO VARIABLES 


This program calculates means, standard deviations, and standard 
errors of the mean for one or two variables, and covariance and 
correlation coefficient for two variables. 


The input is a set of data points: 
1(X1,¥1),(X2,Y2), ... (XVi), »-- (XasVn) } 
in which the y;’s may be omitted if only one variable is being used. 
The statistics computed are: 






x= > X; = mean of x 
i= 1 
‘. 1 n 
y=, > yi = mean ofy 

i=] 

S X;? 7 nx? | H : 

3. — eS ee = Std. dev. of X by n—1 method , @ ben 6 (v 6 oh Bn aS it 
Sx = = std. dev. of x by n method 


= std. dev. of y by n—1 method 

















nA" 
> y? — ny? 
>> = a and = std. dev. of y byn method 
Sx S antan Oe 
> = = Std. error of the mean of x byn—1 method weays or B20 2 
Vn | 
of = Eo .. Std. error of the mean of x by n method 
Vn 
S; = . ae std. error of the mean of y by n—1 method 
Vn | 
oF = ae Std. error of the mean of y by n method 
Vn 
a. = a2 XY; B x x 2 y = covariance of x and y 
‘eh a by n—1 method 
Sky = =| 2 GY; 4 S X; S yi | = covariance of x and y by n method 
i=] i=1 i=1 
— =e 


‘xy = SS = S’ ae correlation coefficient of x and y 


AY 


311-01 


NOTES: 1. The number of data points must be > 1. 

2. This program may be used with the Linear Regression 
program (ST1-08) to determine additional statistics. Data 
points entered on either program may be used for the other. 
Do not initialize if data from the other program is being used. 


Example: Calculate x, y, S,, S, Sy, $5, Sz S& Ss SF, Sig Say, May for the 
data in the following table, then change the x and y values fori = 5 
from (23, 17) to (5, 19) and find the new r,,, and S‘,,. 





OOAON Oat WDM 


ST1-01 


Enter Press Display Comment 
CMs | Initialize 
31 [ A] i XI 
6 1, va 
20.5 [A] 2. X2 
19 Z. y2 
14 [A] 3, X3 
14 3. y3 
etc. 
He 24.38888889 compute x 
~ 2nd) HE =: 16.33333333 compute y 
10.07403152 compute S, 
9.497888004 compute S', 
MM 8 4.981214711 compute S, 
4.696334268 compute S‘, 
[ D | 3.358010507 compute S; 
3.165962668 compute St 
ei 1.660404904 compute Sz 
1.565444756 compute Si 
LE 10.85416667 compute §,,, 
\ 9.648148148 compute Si, 
i EM 2163006968 ’ compute r,, 
23 ue 8. delete 23 (23, 17) 
17 8. delete 17 (23, 17) 
5 [ A] 9. insert 5 (5, 19) 
19 9. insert 19 (5, 19) 
ae 0825567005 compute new r,, 
LE} 5.006944445 compute new S,, 
4.450617284 compute new S, 


ST1-01 
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Sxy, Sxy’ 





USER INSTRUCTIONS 


X) 
‘a 
variable is used X2 
Y2 
Xi 
yi 








TEP 





. 


omit the y steps 
To delete (x;, y;) 

5 Compute outputs 
All outputs except 
those obtained by 


be obtained in 
any order. Outputs 


2) 
= Go} Pm] — 


EEEBEEEEBEEEEBEEEEBREEm 





obtained by pressing 
obtained immediately 


following the 


HAIG Iinin il wm |A Im ia ia t!« | >< | no} po] —a| = 
4 ne E. ‘“|~ |= we ia “an~ | “%~ |x 

< *“< - ; 

4 


corresponding 


keypush shown. 


Bm et ee ee 
EEE EEE 
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46 *LBL 





000 
001 
002 
003 
004 
005 
006 
007 
008 
009 
010 
011 
012 
013 
014 
015 
016 
017 
018 
019 
020 
021 
022 
023 
024 





























“Denotes 2nd function key 


Ca Ce 







11 
42 
00 
01 
ag 
00 
03 
40 
a4 
00 
06 
46 
34 
43 
00 
00 
30 
42 
00 
08 
01 
44 
00 
00 
25 


A 
STO 


PROGRAM LISTING 


56 * 
46 
16 
95 
94 
42 
00 
01 
44 
00 
03 
40 
94 
44 
00 
06 
02 
94 
44 
00 
00 
51 
34 
56 
42 
00 
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at 
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PROGRAM LISTING 
































112 00 O 
113 04 4 
114 41 GTO 
112 So COS 
116 46 *LBL 
WW? SG 
118 43 RCL 
119 00 O 
120 06 6 
127 3. = 
122 53'< 
123 12 6 
124 46 *LBL 
125 24 CE 
126 40 *x? 
127 65 X 
128 43 RCL 
129 00 O 
130 00 O 
131 54 ) 
tse ou * 
133 55 + 
134 43 RCL 
135 00 O 
136 08 8 
137 95 = 
138 56 *rtn 
139 


“Denotes 2nd function key 
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PERMUTATIONS AND COMBINATIONS 


For givennandrwheren =r = 0, the following are calculated. 


Permutations: P(") = a = n-(—1)>... ft) 
eS ON 
~ rl(n—r)! 





Combinations: C(") 


NOTES: 1. After an n,r pair has been entered, either or both may be 
changed for a new calculation. 
2. P(t) = CG) = 1 
3. Error indications 
a. Forr > n, the display flashes 9.999999999 99. 
b. The display will flash 9.999999999 99 in case of overflow in 
the calculation. 
c. For negative entries, the absolute values will be used. 
d. For non-integer entries, the rounded integer value will be 
used. 


Example: Find the number of permutations of 60 items taken 5 at atime 
and the number of combinations of 120 items taken 5 at atime. 


Enter Press Display Comment 
LE Initialize 
60 LA] 60. n 
3) 5. r 
655381440. P(e) 
120 [A ] 120. new n 
[ D | 190578024. (2) 
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TEXAS INSTRUMENTS 


cO-LLS SN.LVNIGINOD 8 SN.LVLNINYAd 2am 
@A PERMUTAT’NS & COMBINAT’NS ST1-02 











Enterr (0 =r =n) 


Calculate P(" 
6) Calculate C(*) 
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PROGRAM LISTING 


084 18 *C’ 
085 42 
086 00 
087 00 
088 43 
089 00 
090 01 
091 42 
092 00 
093 04 
094 75 
095 43 
00 
02 
95 
22 
80 * 
78 
42 
00 
05 
43 
00 
01 
79 
07 
00 


0° Counter OS n— Ff 


—~ answer 
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PROGRAM LISTING 





“Denotes 2nd function key 
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MEANS AND MOMENTS 


For agiven set of input data, {xi, x2,...x,}, with associated frequencies, 
(fi, fo,...f,} (for grouped data), the following means, moments, and 
skewness and kurtosis of distribution are calculated. If f, = 1for alli, the 
calculations are for ungrouped data, otherwise for grouped data. 


Means 
Arithmetic =A =, fx; rarmonte =H =WN/ = 


n \ TIN T/t 
Geometric =G = (11 x") Generalized = M(t) = (a2 f()') 





1 Ni- 
Moments 
he ey ST 

m=A= * FoaPldan m4 = i 2 F:0 A)4 
ee S f(x, — A)? Kurtosis of distribution =—~“ 

N ak i (m2)? 

' S f(x, — A) Sk f distribution = —™3, 
m3 = yj *, (Xx; ) ewness of distribution = (m2)32 
N= S f; n = Number of different groups 


NOTES: 1. The initialization routine assumes a frequency of 1 andt = 1. 

2. The value of t may be changed before any x;’s are entered, if 
changed after any x;, M(t) will be meaningless. 

3. For grouped data if f; = f;_,, it need not be entered. 

4. A negative frequency will cause the next x entered to be 
deleted |f| times. 

5. A negative x will cause the geometric and generalized 
means to be invalid. 


Error Indications: 
A frequency or t of 0 will cause the display to flash 9.999999999 99 and 
the old value to be retained. 


Example: The ungrouped (f, = 1) data {x; = 1.1, 2.3, 4.7, 3.5, 1.8, 2.9} 
yields the following: A= 2.716666667, G = 2.451843924, H = 2.187595899, 
m(1) = 2.716666667; m2 = 1.368055556, m3 = .5367592591, 

ma = 3.903019215; kurtosis of distribution = 2.085418497, skewness 

of distribution = .3354470243 


References: Handbook of Mathematical Tables and Formulas, R. S. 
Burington, McGraw-Hill, Fourth Edition, 1965, p. 164. International 
Dictionary of Applied Mathematics, W. F. Frieberger, D. Van Nostrand, 
1960. 
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L-€0-1LS (L)SLNSWOW 8 SNVSN 2am | 
|| @Am MEANS & MOMENTS (1 ST1-03-1 
SCOT TSE RE ee Ge Se 


c-€0-LLS (Z)SLNSWOW 8 SNVAW 2am | 
|  @Am MEANS & MOMENTS (2 ST1-03-2 
Ae. CURE 5 Oras ee an: Be et Gee ae 
kurt 
USER INSTRUCTIONS 


Enter Program ST1-03-1 





skew 


NO 


3 Enter t (for gen. mean) 


a 


Ww 
=| 
Uv 


Enter data 
a. for grouped data 


Cnt. BT, He he i, #0 


t4#0 
Xi 
Xj 


corrected Sf 


d. kurtosis 


f. means, see step 6 


—™~] 


~~ 


urtosis 
Skewness 


Pelle] | elelele) fale! | | ele) | ee 
LE Ef elmer 
BERR RRRRERRRRRRR EERE 


seatl & 
ee) 
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PROGRAM LISTING 





46 *LBL 
11 A 


028 
029 
030 
031 
032 
033 
034 
035 
036 
037 
038 
039 
040 
041 
042 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 
053 
054 





























055 65 X 


084 00 0 
085 05 5 
OS6 55 + 
087 
088 
089 
090 
091 
092 
093 
094 
095 
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PROGRAM LISTING 


© 


140 11 A 46 iy | 
141 01 1 *C: ge *1/x 
142 94 43 RCL 56 *rtn 
143 22 00 0. 00 
144 49 01 90 
145 00 55 - 00 
146 02 43 00 
147 43 00 00 
148 00 07 00 
149 01 95 = 00 
150 56 56 00 
151 46 46 00 
152 16 19 00 
153 43 43 00 
154 00 00 
155 08 00 
156 55 + 00 
157 43 00 
158 00 00 
159 01 00 
160 95 00 
161 35 V 00 
162 43 00 
163 00 00 
164 03 00 
165 95 = 00 
166 56 00 
167 46 00 





“Denotes 2nd function key 


S1T1-03- 2 


PROGRAM LISTING 





















































































000 46 *LBL 056 55 + 084 09 9 
001 11 A 085 55 + 
002 43 RCL 086 43 RCL 
003 01 1 087 00 0 
004 02 2 088 01 1 
005 40 *x? 089 75 — 
006 94 +/-— 090 03 3 
007 85 + 091 65 X 
008 53 ( 092 43 RCL 
009 43 RCL 093 01 1 
010 00 O 094 02 2 
011 09 9 095 40 *x? 
O12 55 + 096 40 *x? 
013 43 RCL O97 95 = 
014 00 O 098 56 *rtn 
O15 01.1 099 46 *LBL 
016 95 = 100 14 D 
017 56 *rtn 101 11 A 
018 46 *LBL 102 40 *x? 
019 12 B 103 42 STO 
43 104 01 1 
105 03 3 
106 13 C 
107 46 *LBL 
108 87 *1’ 
109 55 + 
110 43 RCL 


111 


“Denotes 2nd function key 
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PROGRAM LISTING 



























112 03 


*Denotes 2nd function key 


3 














































140 35 Vy 


113 95 = 141 43 RCL 
114 56 *rtn 142 00 O 
115 46 *LBL 143 01 1 
1146 15 E 144 95 = 
17At & 145 56 *rtn 
118 65 X 146 46 *LBL 
119 40 *x2 147 18°C’ 
120 95 = 148 43 RCL 
121 30 *Vx 149 00 0 
122 95 = 150 01 1 
123 42 STO 151 55 + 
124 01 1 152 43 RCL 
125 03 3 153 00 0 
126 12 B 154 07 7 
127 41 GTO 155 95 = 
128 87. *1’ 156 56 *rtn 
129 46 *LBL 157 46 *LBL 
130 16 *A’ 158 19 *D’ 
131 43 RCL 159 43 RCL 
132 01 1 160 00 0 
133 02 2 161 08 8 
134 56 *rtn 162 55 = 
135 46 *LBL 163 43 RCL 
136 17 *B’ 164 00 0 
137 43 RCL 165 01 1 
138 00 0 166 95 = 
139 167 35 Vy 
































S1T1-04 


RANDOM NUMBER GENERATOR 


This program generates uniformly distributed random numbers 
(0 = U; = 1) by use of the following formula: 


Ui,, = 77U,;(Mod 107°) 
Where: U,, = 1 


Normally distributed random numbers are also generated, using the 
direct method: 


A; = ( — 2th uy" COS (Q7rU, . t) 
Ais, =(—2In U;)'’? sin (27U; , 1) 


Where: U; and U,, , are generated from the above method. 
The final result is: 


Y,=oA,+m 


which yields normally distributed random numbers with mean m and 
standard deviation o. The user may input the standard deviation, mean, 
and U,, or the default values are: U; = 0.55, 7 = 1,andm=0O. 


NOTES: 1. Fornormal distribution, the calculator must be in radian mode. 

2. The random numbers return as fixed point 5. 

3. The initial seed is 1.3597878 x 10'°, and is entered when [ A |] 
is pressed. 

4. The [e—] and Ww are of use only with the printer. [ E | 
initiates continuous printing of U;, and Mw continuous 
printing of Y;. Printing is stopped by pressing and then 
a . 


Example: To generate uniformly distributed random numbers the 
procedure is: 


Enter Press Display Comment 
LA | 0.00000 Initialize 
0.48405 U; 
0.46815 U2 


Each pressing of will display a new U,. With a printer, pressing [E | 
once instead of will produce a list of U; which is stopped by 
pressing and then ' 


To generate normally distributed random numbers the procedure is: 


Enter Press Display Comment 
LA | 0.00000 Initialize 
2 2.00000 mean 
1 Ae 1.00000 o 
[ D | 0.91202 Y 
[ D | 2.10940 Y2 


ST1-04 


Each pressing of LD _] will display a new Y;. With a printer, pressing 

Mal once instead of [_D ] will produce a list of Y; which is stopped by 
pressing and then _If the mean and o are not entered, the 
default values are m = Oando = 1. 


References: Theory and Problems of Numerical Analysis, Schaum's 
Outline Series, McGraw-Hill Book Company. Handbook of Mathematical 
Functions, U.S. Department of Commerce, National Bureau of 
Standards, 1964. 
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USER INSTRUCTIONS 


a 

Mi | exerProgram aanaay [st dC 

Te [sideomswicn | | iY 

Ds 2k ILS AA, LO 

8 a 

ra [compute it —S—«d a 

[each pusnproduees 

es eT 

| entermean mean ‘|e 
ee 


















ia 
i 
a 
a ee 
—}+—_— 
Ewnareliet 
| Eachpush produces | hee 
Wes See a 
8 | PrintertistofU; | ii 
9 | Printertistofry, vom 
10 stopstepsoro | 
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000 46 *LBL 
















001 
002 
003 
004 
005 
006 
007 
008 
009 
010 
011 
012 
013 
014 
015 
016 
017 
018 
019 
020 
021 
022 
023 
024 
025 
026 
027 


“Denotes 2nd function key 


11 A 
01 
93 
03 
05 
09 
07 
08 
07 
08 
of EE 
O1 1 

00 0 


—. 


CON ON O OW 


22 INV 


oe EE 
42 
00 0 
00 0 


22 INV 
50 “st fig 


O1 1 
93 « 
US'S 
05 5 


42 STO 


00 0 
06 6 


PROGRAM LISTING 


028 01 1 
STO 


056 
057 
058 
059 
060 
061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
071 
072 
073 
































AY 
00 
09 
95 
ex 
o2 
és 
93 
O05 
o/ 
00 
95 
92 
22 
92 
sf 
05 
65 
01 
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PROGRAM LISTING 












































































































































112 00 0 140 42 STO 196 00 0 
113 04 4 141 01 1 197 01 1 
114 93 e 142 00 0 198 16 *A’ 
115 055 143 33 cos 199 56 *rtn 
116 05 5 144 42 STO 200 46 *LBL 
117 42 STO 145 00 0 201 15 E 
118 00 0 146 07 7 202 12 B 
119 06 6 147 43 RCL 203 98 “prt 
120 43 RCL 148 00 0 204 41 GTO 
121 00 0 149 06 6 205 15 £ 
122 04 4 150 23 Inx 206 46 *LBL 
123 56 *rtn 151 65 X 207 10 *E’ 
124 46 *LBL 152 02 2 208 14 D 
125 14D 153 94 +/- 209 98 “prt 
126 60 “if flg 154 95 = 210 41 GTO 
127 01 1 155 30 *Vx 211 10 *E’ 
128 68 *8’ 156 42 STO 212 00 

129 50 *stflg | 157000 213 00 

130 01 1 158 01 1 214 00 

131 12B 159 16 *A’ 215 00 

132 59 *z 160 56 *rtn 216 00 

133 65 X 161 46 *LBL 217 00 

134 43 RCL 162 16 *A’ 218 00 

135 00 0 163 65 X 219 00 

136 08 8 164 43 RCL 220 00 

137 65 X 165 00 0 221 

138 02 2 166 07 7 





139 95 = 167 65 X 





“Denotes 2nd function key 
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ONE-WAY ANALYSIS OF VARIANCE (ANOVA) 


One way analysis of variance is a technique for testing the differences 
between the population means of k treatment groups, where each group 
i(i = 1, 2,...,k) consists of n; observations x,,(j = 1, 2,...,n;). The 

different groups need not have the same number of observations. 


Formulas used: 
nj 


Sum of observations in groupi= >» Xj; 
j=1 


k nj F 
cS Pe > x: 
Total Sum of Squares = SS= > & x2 —-———— 
i=] j=1 . n, 
nj ach kK Nj 
k (> xi} (= 2 x.) 
Treatment Sum of Squares = TSS = 2 — — MELi=1 
- | Sn, 


Error Sum of Squares = ESS = SS — TSS 
The value of the F statistic is: 


. Ore x TSS 
~ DFi x ESS 


Where: DFi = Treatment degrees of freedom = k — 1 


F 


k 
DF2 = Error degrees of freedom = > n,—k 


i=! 


X; Oi;" 

1) 25 20°12 22 ie 18.66666667 }#18.55555556 
2130 45 25 46 19 33 37 33.57142857 85.10204082 
Sin® 19 14° 764-7 48 21 12.14285714  43.55102041 


Example: 





F = 14.11952 
DFi =2 
UF2= 17 
SS = 2692.8 


TSS = 1680.895238 
ESS = 1011.904762 
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S0-11S JONVIYVA JO TWNV AVM-SNO 2am I | 
[| #Am ONE-WAY ANAL OF VARIANCE ST1-05 
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USER INSTRUCTIONS 


Puen aban | ren pr) 





TEP 


ho 


wm 
Cw —_ 





ei | 
i ea) 
eer) 
ee 
! ort aes 
| (must bedonebefore | 
psteps) tm 
5 | Calculator? | 
(Step Smustbe | 
[performed foreachi) | 

rel as: 

bo, OU Eo 

Tete 


Display o-;? 


eae 
8 | Calculate F,DF1,DF2 re Gal. k 





“TT 


eles] | fe) | eel | BB 
BREE ERR RRe 
A 





ee | Display Treatment DF Bo ety 4 DFi 
Display Error DF eae We DF2 
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PROGRAM LISTING 


: 028 00 0 084 00 0 
029 03 3 085 01 1 
030 40 *x? 086 95 = 
031 94 +/- 087 42 
032 44 SUM 088 01 
033 00 0 089 09 
034 04 4 090 43 
035 01 1 091 00 
036 94 +/-— 092 04 
037 44 SUM 093 44 
038 00 0 094 00 
039 01 1 095 06 
040 56 *rtn 00 
041 46 *LBL 42 
042 13 C 00 
043 01 1 01 
044 44 SUM 42 
045 00 O 00 
046 00 0 03 
047 43 RCL 42 
048 00 0 00 
049 01 1 04 
050 44 SUM 43 
051 00 O 01 
052 02 2 08 
053 43 RCL 56 
054 00 0 
055 03.3 
































7 X(2Xij)2/n; 
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PROGRAM LISTING 


"Denotes 2nd function key 


Zo 


ST1-06 


t STATISTIC EVALUATION 


(Paired Observations) 


This program evaluates the t statistic 


5 


_AVn 
t ee. 
with n—1 degrees of freedom, to test the difference between the means 


(unknown) of two normally distributed populations, using paired values 
of data from the two populations. 





Given the set of n pairs of x,y;, let 


iby = Xp Ys; 
yeh n 
A = = ba A; 
$ a2—1(3 4)" 
Sa — | i= ; n i=1 ) | 
iia 
X; 6 4 7 4 3 
Example: y, 5 3 5 4 3 
t = 1.870828693 
DF = 4. 
A= 14 
S, = 1.673320053 
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90-LLS G3auIVd) WAS LVLS-} =a | 
|| @Am t-STAT EVAL (PAIRED ST1-06 
ce ae ee a ee es ee 





USER INSTRUCTIONS 















Pt | EnterProgram(aands) [| TTT 

9 SER (eee WS) ee a 
ie. aaa 

To delete data eae a 

Pentex |r) | 

[entry || 

| 6 | ComputeDF | CC 
| Computes: | | En) 


2/ 
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PROGRAM LISTING 


AY 
000 56 *rtn 43 









































































001 46 *LBL 057 00 0 
002 11 A 058 03 3 
003 42 STO 059 56 *rtn 
004 00 0 060 46 *LBL 
005 01 1 061 13 C 
006 56 *rtn 062 43 RCL 
007 46 *LBL 063 00 0 
008 12 B 064 02 2 
009 95 = 065 55 + 
010 94 +/- 066 43 RCL 
011 85 + 067 00 0 
012 43 RCL 068 03 3 
013 00 O 069 95 = 
014 01 1 070 42 STO 
O15 95 = 071 00 O 
016 44 SUM 072 04 4 
017 00 0 073 40 *x? 
018 02 2 074 65 X 
019 40 *x? 075 43 RCL 
020 44 SUM 076 00 0 
021 00 0 O77 03 S 


08 











94 
053 44 SUM 
054 00 0 
055 03 3 
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PROGRAM LISTING 


140 15 E 
141-25 GLA 
142 47 *CMs 
143 86 *rset 
144 00 

145 00 

146 00 

147 00 

148 00 

149 00 

150 00 

151 00 

152 00 

153 00 

154 00 

155 00 

156 00 

157 00 

158 00 

159 00 

160 00 

161 00 

162 00 

163 00 

164 00 

165 00 

166 00 

167 00 





“Denotes 2nd function key 
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t STATISTIC EVALUATION 
(TWO SAMPLE TEST) 


This program evaluates the t statistic 


247s 
(+ ag 1)" (x Se nx? nibs Ly; eo ny 
fh, N, Ph, =a 


with n, +n, — 2 degrees of freedom, to test the significance of the 
difference (A) between the means (unknown) of two independent random 
samples from two normally distributed populations having the same 
unknown variance. 


= 


Given the sample: 


i EN, MS, on KI 
¥i = [Y1, ¥2: .... Vay] 
then gus ss x, 
Nx i=] 
= Ne Ny 
y n,, Pe yj 


Example: x; : 71, 82, 105, 110, 124, 111, 126, 129 
y; >o7, 75,66, 93, 98, 102, 115, 117, 121,180 


For A=0, t= .7562582916 
DF = 16. 
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LO-LLS “(A1d VS OML) 1WAS LVLS-} =a | 
| @A t-STAT EVAL (TWO SAMPLE ST1-07 
ae ee eee oa eae ee eee 











ae 
a ee eat 
mem b!  e ornt 
eC. Miia bs sig 
2 ae eT 
25 a 
Eo ET 
, i eae 
ee Lt 
idbeoks nl © RIN ho: 
‘he Bh 1 eee 
a ee eee 
ee i ea 
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PROGRAM LISTING 


56 *rtn 084 00 O 
46 *LBL 085 08 8 
11 A O86 40 *x? 
087 85 + 
088 53 ( 
089 43 RCL 
090 00 
091 05 
092 75 
093 43 
094 00 
095 06 
65 
43 
00 
09 
40 
95 
os 


43 
01 
00 
95 = 
30 
42 
01 





311-QOF 


PROGRAM LISTING 


112 00 O 
113 03 3 
114 20 *1/x 
113 85 
116 43 
117 OO 
118 06 
118 20 = 
120 95 
121 30 
lez 6S 
123 43 
124 01 
120 O01 
12635 = 
127 20 
128 65 
129 53 
130 43 
131 

132 

133 

134 

135 

136 

137 

138 

139 





“Denotes 2nd function key 
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| | LINEAR REGRESSION 
The linear least-squares fit of input data points (x,y) is calculated using 
the following: 
XD xy 
Slope = m= i 
= Pe 
n 
ZY = max 


Intercept = b = — = 


The y’ for entered x and x’ for entered y are calculated as follows: 
y’=mx+b 
x’ =(y—b)/m 


The coefficient of determination is calculated as: 
a. 


NOTES: 1. The calculator must be initialized by pressing [IS before 
each new set of data is entered. 

2. The number of data points entered must be greater than 1 
before any of the calculations are attempted. 

3. This program may be used with the Basic Statistics program 
(ST1-01) to determine additional statistics. Data points 
entered on either program may be used for the other. Do not 
initialize if data from the other program is being used. 


Example: 

Kizo.4 3:7 218° 20.1 20.7 
Yi 92 B68: 114 34° 03 
r? = .9256678121 

m= 1.22906793 

b =—15.40505529 


Fory = 12, x’ = 22.29742931 
For x = 22.3, y’ = 12.00315956 
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x—>y | «| delete [| 
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Enter Prosram (Rand) [ 
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PROGRAM LISTING 


































46 *LBL 
19 *D’ 
95 = 
94 +/— 
42 STO 
00 O 

02 2 

44 SUM 
00 0 

04 4 

40 *x? 
94 +/— 
44 SUM 
00 0 

O7 7 

43 RCL 
00 0 

02 2 

65 X 
43 RCL 
00 0 

O01 1 

95 = 
44 SUM 
00 0 

05 5 

43 RCL 
00 0 


aN -08 


\ 


PROGRAM LISTING 
bia Take LT 


140 77 *4’ 
141 42 STO 
142 01 1 
143 02 2 
144 46 *LBL 
145 15 
146 14 
147 65 
148 43 
149 00 
150 03 
151 
152 8&5 
153 43 
154 00 
155 04 
156 95 = 
S67 Sa 
158 43 
159 00 
160 O00 
161 95 = 
162 56 
163 46 

O01 1 164 16 

O1 1 165 51 

56 *rtn 166 77 

46 *LBL., | 167 94 


“Denotes 2nd function key 
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S1T1-09 
f 


| y POWER CURVE FIT 
dive a set of data points: 
{XiYi); )= 1, 2 é seis n} 
Where: x; > 0 and y, > 0, this program will fit a power curve: 
y = ax? (a > 0) 


The problem may be interpreted as a linear regression problem by using 
the equation in the form: 


Iny=bInx-+Ina 


The computed statistics are: 
1. Regression coefficients 


—_— E Iny, b>In A 
n n | 
E (in x¢in yi) — EE TO yy) 
b ee ee a ee ee eee Se oe 
— (& In x,)? 
X(In_ x;)? eo 


NOTE: nis a positive integer + 1. 
2. Correlation coefficient 


S (In x,)(In y,) Gin ee ; 


b in. xjF— caer b (In y;)? — ane 


r2 = 


3. Estimated value x’ for given y 


b 
x'= Vy/a 
4. Estimated value y’ for given x 
y’ = ax? 


11 13.0 15: "B85 s920 
oO: TABOR O89. 188: 2.4 ore Sa 
a = .2901902009 
= ./759169351 
r? = 0.921727263 
forx=9 , y’=1.596230726 
for y = 2.8, x’ = 18.56904361 





















Example: 
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y ‘ 

4 
me 
7 & 
2 e 
a 
& 

1 & 
0 x 

0 5 10 15 20 





USER INSTRUCTIONS 


STEP PROCEDURE ENTER PRESS DISPLAY 










4.] EntorPrepram (aand@)” |). Pee ee N | eee 
initialize ae ee 0° ae 
3 | Enterdata CD 
|| repeatforal) sf 


gs Delete data Pe ae ae aa 

dick Tae 

Pf bemery fy 

[| eveticints TN 

6 | Caleuateconeation | 

Sk cee 

C7 | Computeestimated | 
i eae 
Pee 


accel 
EF 
age 


\ 
\ 
‘ 











value for 
x’ given y 
y’ given x ——= 


tele | | (elelslele 
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PROGRAM LISTING 


000 46 *LBL 056 46 *LBL 

























































001 11 A O57 16 *A’ 
002 47 *CMs 058 23: inx 
003 25 CLR 059 42 STO 
004 56 *rtn 060 00 0 
005 46 *LBL 061 07 7 
006 12 B 062 22 INV 
007 23 Inx 063 44 SUM 
008 42 STO 064 00 0 
009 00 0 065 03 3 
010 06 6 O66 40 *x? 
011 22 INV 067 22 INV 
012 23 Inx 068 44 SUM 
013 56 *rtn 069 00 0 
014 46 *LBL O70 05-5 
015 13 € 071 43 RCL 
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PROGRAM LISTING 






























































“Denotes 2nd function key 































4] 

































20 *1/x 140 08 8 168 46 *LBL 196 01 1 
113 42 STO 141 94 +/— 169 17 *B’ 197 00 0 
114 01 1 142 65 X 170 55 = 198 65 X 
115 00 0 143 43 RCL 171 43 RCL 199 53 ( 
116 65 X 144 00 0 172 00 0 200 43 RCL 
117 53 ( 145 01 1 173 09 9 201 00 0 
118 43 RCL 146 85 + 174 95 = 202 05 5 
119 00 0 147 43 RCL 175 35 Vy 203 75 — 
120 04 4 148 00 0 176 43 RCL 204 43 RCL 
121 75 — 149 03 3 177 00 0 205 00 0 
122 43 RCL 150 95 = 178 08 8 206 03 3 
123 00 0 151 55 + 179 95 = 207 40 *x? 
124 01 1 152 43 RCL 180 56 *rtn 208 55 + 
125 65 X 153 00 0 181 46 *LBL 209 43 RCL 

43 154 00 0 182 18°*C’ 210 00 0 
155 95 = 183 45 y* 211 00 0 
128 03 3 156 22 INV 184 43 RCL 212 54 ) 
129 55 + 157 23 Inx 185 00 0 246 80 *isx 
43 158 42 STO 186 08 8 214 65 X 
159 00 O 187 65 X 215 43 RCL 
160 09 9 188 43 RCL 216 01 1 
161 56 *rtn 189 00 0 217 01 1 
162 46 *LBL 190 09 9 218 40 *x2 
163 19 *D’ 191 95 = 219 95 = 
164 43 RCL 192 56 *rtn 220 42 STO 
165 00 0 193 46 *LBL 221 01 1 
166 08 8 194 15 E 222 02 2 
167 195 43 223 56 “Tih 
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| EXPONENTIAL CURVE FIT 
iven a Set of data points: 
1(Xi; Yi), § = 1, 2,..., n} y, >0 


This program finds the least-squares fit for an exponential function of 
the form: 


y = ae’* (a > 0) 
The problem is computed according to the linear equation 
Iny = Ina + bx 


Statistical outputs: 
1. Coefficients a, b 


a = exp pe bey 


) = ————— eee 
| sy2 — (Xi)? » 
3 n 


Note: nis positive integer + 1. 
2. Correlation coefficient 


} -|2 
ExiIny Get amen 


fa 43 NS 
3x7 | S(Iny,) _ (lny)) | 
n n 
3. Estimated value for x’ given y | 


4. Estimated value for y’ given x 


y’ — aes 


Example: X; 1.42 1.94 2.63 3.11 a 4.25 





yi .63 9 1.42 1.9 2.43 2.89 


a= .3131935086 
b= .5477927735 
r?= .9810626986 


Forx=1.6, y' = .7524162355 
Fory=2.2, x’ =3.558629277 
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0 1 2 3 4 2 






{| _ @Am™ EXPONENTIAL CURVE FIT ST1-10 
| delete | y—>x | x—y [| bb | 
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PROGRAM LISTING 






































“Denotes 2nd function key 












































































000 46 *LBL 028 00 0 056 56 *rtn 084 44 SUM 
001 10 *E’ 029 05 5 057 46 *LBL 085 00 0 
002 43 RCL 030 43 RCL 058 16 *A’ 086 02 2 
003 00 0 031 00 0 059 23 Inx 087 43 RCL 
004 00 0 032 06 6 060 42 STO 088 00 0 
005 56 *rtn 033 44 SUM 061 00 0 089 O07 7 
006 46 “LBL 034 00 0 090 22 INV 
007 11 A 035 01 1 091 44 SUM 
008 47 *CMs 036 49 *PROD 092 00 0 
009 25 CLR 037 00 0 093 04 4 
010 56 *rtn 038 07 7 094 01 1 
011 46 *LBL 039 40 *x? 095 94 +/-— 
Ole te 6 040 44 SUM 096 41 GTO 
013 42 STO 041 00 0 O37 87 *1’ 
014 00 0 042 02 2 098 46 *LBL 
015 06 6 043 43 RCL 099 14 D 
016 56 *rtn 044 00 0 100 43 RCL 
017 46 *LBL 045 07 7 101 00 O 
018 13 C 046 44 SUM 102 02 2 
019 23 Inx 047 00 0 We 75 > 
020 42 STO 048 04 4 104 43 RCL 
021 00 0 049 01 1 

022 O07 7 050 46 *LBL 

023 44 SUM O51 67: *t’ 

024 00 0 052 44 SUM 

025 03 3 053 00 0 

026 40 *x? 054 00 0 95 = 
027 44 SUM 055’ 70 TE" V11 eer 1 /x 
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PROGRAM LISTING 



































196 00 0 
197 65 X 
198 53 ( 
199 43 RCL 
200 00 0 
201 05 5 
202 75 — 
203 43 RCL 
204 00 0 
205 03 3 
206 40 *x? 
207 55 + 
208 10 *E’ 
209 54 ) 
210 20 *1/x 
211 65 X 
212 43 RCL 
213 01 1 
214 01 1 
215 40 *x? 
216 95 = 
217 42 STO 
218 01 1 
219 02 2 
220 56 *rtn 
221 00 

222 00 
223 00 


“Denotes 2nd function key 
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LOGARITHMIC CURVE FIT 
Given a Set of data points: 
{(X;, Yi), |= 1, 2,...,n} 
Where: x,;>0, this program fits a logarithmic curve: 
y=a-+ binx 


Statistical outputs: 
1. Regression coefficients 


ZV, — O21nXx, 


a = 





3. Estimated value for x’ given y 
x’ = exp b b | 


4. Estimated value for y’ given x 





y’=a+blinx 
NOTE: nis positive integer # 1. 


Example: 





yi 9.3 10.8 Wo 28.2 36.9 44.4 
a = —13.49953513 


b= 18.91052316 
r? = 8418448418 
For x = 21, y’ = 44.07397695 


For y = 40, x'= 16.93001448 
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LL-LLS Lid SAUND DINHLINVDO1 =a || 
||  @A® LOGARITHMIC CURVE FIT ST1-11 
vane eee 8! be Ee 


delete —> x’ x— y’ 





USER INSTRUCTIONS 
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PROGRAM LISTING 


0 084 00 O 
085 07 7 
086 22 INV 
087 44 
088 00 
089 04 
090 O01 
091 94 
O92 41 
093 87 
094 46 
095 14 

43 
00 
02 






































PROGRAM LISTING 





























































































































































112 00 0 140 43 168 65 X 
113 65 X 141 00 0 169 53 ( 
114 53 ( 142 01 1 170 43 RCL 
115 43 RCL 143 85 + 171 00 0 
116 00 0 144 43 RCL 172 05 5 
117 04 4 145 00 0 173 75 — 
118 75 — 146 03 3 174 43 RCL 
119 43 RCL 147 95 = 175 00 0 
120 00 0 148 55 + 176 03 3 
121 01 1 149 43 RCL 177 40 *x? 
122 65 X 150 00 0 178 55 + 
123 43 RCL 151 00 0 179 43 RCL 
124 00 0 152 95 = 180 00 0 
125 03 3 153 42 STO 181 00 0 
126 55 + 154 00 0 182 54 ) 
127 43 RCL 155 09 9 183 20 *1/x 
128 00 0 156 56 *rtn 184 65 X 
129 00 0 157 46 *LBL 185 43 RCL 
130 54 ) 158 19 *D’ 186 01 1 
131 42 STO 159 43 RCL 187 01 1 
132 01 1 160 00 O 188 40 *x? 
143 01 1 161 08 8 189 95 = 
134 95 = 162 56 *rtn 1190 42 STO 
135 42 STO 163 46 *LBL 191 01 1 
136 00 0 164 15 E 192 02 2 
137 08 8 165 43 RCL 193 56 *rtn 
138 94 +/— 166 01 1 194 46 *LBL 
139 65 X 167 00 0 195 





*Denotes 2nd function key 
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MULTIPLE LINEAR REGRESSION 


Using the least-squares method this program (2 cards) derives a linear 
equation of the form: 


Z=ao+ aix + azy 
for a given set of data points: 
1(X), Yi. Zi), |= 1, 2,...,N} 
According to the normal equations 
2Z; = aon + aidx; + a2dy; 
2X;Zj = AoLX; + aiLx? + a2dx,y; 
2ZYiZ; = Aody; + aidxiy; + a2dy? 
the regression coefficients ao, a1, a2 may be computed. 
+ Ni — No 
[N&x? — (2x;)?] [n2y? — (Zyi)?] — [nexiy; — (2x;)(Zy;) ]? 
Where: Ni = [n&x? — (2x;)?] [nZy,z, — (Ly;)(=z;) ] 
N2= [n2xiy; — (2x;)(Zy;)] [n2x,z; — (2x;)(2z;) ] 
_ [nxiz, — (2x;)(2z;) ] — a?[n=xy; — (2x;)(Zy;) ] 
Nike — (2x;)* 
22) > A2LY; => ALEX, 
n 


a2 


al 


ao = 


NOTE: nis positive integer + 1 


Card 1 is used when entering data points. If an error is made while 
entering coordinates, the incorrect data may be deleted by entering x, 
and y; in the normal manner. Then, z, should be entered in the display 
and [_D | pressed to delete x, y, and z. When all data points have been 
entered, any sum may be recalled by pressing and the corresponding 
data register number (refer to User Instructions). 


Card 2 must be loaded before calculating the regression coefficients. 
After the coefficients are determined, the estimated value for z,(z’) may 
be computed for a given x; and y;. | 


Reference: Computer Methods for Science and Engineering, Robert 
LaFara, Hayden Book Co., 1973 
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1.3 
1.24 
4.4 


2.6 
29 


Example: 













1.88 








1. Regression line: 
Z = .6462676006 + .4986757882x + 1.242837142y 
2. Regression coefficients: 
ao = .6462676006 


ai = .4986757882 

a2 = 1.242837142 

For: x= 1, y= 2.5, z' = 4.252036243 
Registers: Roo n=4 Ros yi? = 14.692 
Ro 2X: — 6.65 Ros &XY; = 14.113 
Ro2 2y; = 7.12 Foy 2%2, = 28.875 
Pos..3,2Z, = 14.65 Ros Xy;Z; = 29.899 


Ros 2X? = 13.8525 
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@A® M’TIPLE LIN REGRESSION (2 ST1-12-2 
RM Se cee Dee ee 
ao ai a2 x y=—>z’ 
USER INSTRUCTIONS 
Ca a a ae 
ee 
eee 
[a ene ata [xa i. 
tran ter 
Te i re 
a 
SS a CR 2] 
a) 
Rk CMA TARA =| ea 
= [esas | [wear] COT Can 
ae be recalled, if desired Yee RCL raj) ae 
a (2 TNs 
a NN (2) Ce Ce 
rr 
ee UCN 
[ronan sree | 
7 ca egesin [Ta = 
eis 
1 Eee ET |. +! a 
Fe i haces EE CS | SR a 
Li | wivetorzghenxandy [yo tee 
Lie | (must be done after step 7) ew 6 5 Obs See 
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PROGRAM LISTING 





46 
15 
47 
20 
96 
46 
006 11 
42 
01 
02 
96 
46 
12 
42 
01 
03 
96 
46 
13 
44 
00 
03 
42 
01 
05 
42 



























*L.BL 
E 
*CMs 
CLR 
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01 
03 
49 

























04 
44 
00 
02 
42 
01 
06 
40 
4a 
00 
05 
43 
01 
02 
49 
01 
05 
49 
01 
06 
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*Denotes 2nd function key 
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084 43 RCL 
085 01 1 
086 05 5 
087 65 X 





000 
001 
002 
003 
004 
005 
006 
007 
008 
009 
010 
011 
012 




















46 *LBL 


11 A 


43 RCL 


00 0 
00 0 
65 X 


43 RCL 


00 O 
04 4 
io. = 


43 RCL 


00 0 
O1 1 


40 *x? 


95 = 


42 STO 


01 


PROGRAM LISTING 





















03 
43 
01 
04 
65 
93 
43 
00 
00 
65 
43 
00 


20 *1/x 
70 *iferr 
97 “list 
65 X 


( 
RCL 









088 
089 
090 
091 
092 
093 
094 
095 
096 
097 
098 
099 
100 






























101 








102 








103 








104 








105 








106 


















107 
108 
109 
110 


93 
43 
00 
00 
65 
43 
00 
07 
¥5 
43 
00 
01 
65 
43 
00 
03 
94 
42 
01 


( 
RCL 










































113 
114 
119 
116 
tal 
118 
119 
120 
121 
122 
123 
124 
129 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 


“Denotes 2nd function key 


112 65 X 


43 
01 
O05 
94 
85 
43 
01 
06 
JS 


RCL 
; 

5 
+/- 
—f 
RCL 


o> 


43 
01 
04 
JO 
42 
01 
00 
65 
43 
00 


PROGRAM LISTING 


196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
Zi 
216 
217 
218 
21g 
220 
221 
222 
223 
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09 9 
95 = 


41 GTO 


14 D 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
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QUADRATIC REGRESSION 


The quadratic equation 
Y = a2X? + aiX + Ao 
with normal equations 


n n n 
~ ¥;= ard) kP* +a & x, + aon 
i=] i=] 


t=] j= i= 


n 
ys XiVj = a2 P x9 + ai 2 X;? + ao 
= i 


1 i=1 i= 


n n n n 
D x;"*¥i— ao > xP +a 2 xy Peo 2 x; 
- i=1 i=1: i= 


is used to approximate functions of a single variable whose true form is 


not known, where 
n = number of (x,y) points (n = 3,4,5,...) 
Also, for a given x, this program will calculate 
y(x) = a2x? + aiX + ao 
Where: a2, ai, and ao are determined from the normal equations. 


NOTES: 1. To delete an x,y pair, enter x on[_ A ]andyon . The 
displayed count, n, will be reduced by 1. If you attempt to 


delete an x,y pair before data has been entered, —1 will flash 


in the display. © 


2. If x; is not entered, it will automatically be assigned the value 


of n when y;, is entered. 


Reference: Computer Methods for Science and Engineering, Robert L. 


LaFara, Hayden Book Co., 1973 
Example: Given the sample: 


X 0 1 2 3 A =, 6 7 8 
y 3.2 2 T.2 a 1 2 Be 4 6.6 


find a2, ai, ao, and y(6.2). 


a2 = .2392857143 
ai = —1.505952381 
ao = 3.2 


Thus, y= .2392857143x?2 — 1.505952381x + 3.2 
and y(6.2) = 3.061238095 
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PROGRAM LISTING 


24 GE 
46 *LBL 
52 &E 
26. * 44x 
OF)-1 
94 +/— 
56 “rtn 
00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 





“Denotes 2nd function key 
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PROGRAM LISTING 





























000 46 *LBL O56 20 “1/x 
001 11 A 057 70 *iferr 
002 43 RCL 058 97 “list 
003 00 0 059 65 X 
004 00 0 060 53 ( 

005 65 X 061 43 RCL 
006 43 RCL 062 O1 1 


007 00 O 
008 02 2 
O09 7/5 = 
010 43 RCL 
011 00 O 
012 O1 1 
013 40 *x? 
014 95 = 
015 42 STO 
016 O1 1 
017 04 4 
018 65 X 
019 53 ( 

43 


















“Denotes 2nd function key 


REGISTERS 


Ls C2 eS CE 
13 
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PROGRAM LISTING 






DISPLAY 


































































































































112 65 X 168 13 C 
113 43 RCL 169 43 RCL 
114 01 1 170 O1 1 
115 05.5 we O44 
116 94 +/— 172 56 *rtn 
117 SS + 173 46 “LBL 
118 43 RCL 174 14 D 
17901... 175 42 STO 
120 06 6 176 O1 1 
121 95 = 177 02 2 
122 55> 178 40 *x? 
123 43 RCL 179 65 X 
124 01 1 180 43 RCL 
125 04 4 181 O01 1 
126 95 = 182 01 1 
127 42 STO 183 85 + 
128 O01 1 184 43 RCL 
129 00 0 185 O1 1 
130 65 X 186 02 2 
131 43 RCL 187 65 X 
132 00 0 188 43 RCL 
133 01 1 189 O01 1 
134 94 +/— 12 190 00 O 
135 85 + 181 66.4 
136 43 RCL 164 01 1 192 43 RCL 
137 00 O 165 00 O 193 00 O 
138 05 5 166 56 *rtn 194 09 9 
138 1S 46 *LBL 195 95 = 





“Denotes 2nd function key 
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HISTOGRAM CONSTRUCTION 


This program constructs a histogram, given the range (Xin tO Xmax) and 
the individual observed values (x;) of a variable. Each data value entered 
is checked against the range, and those meeting the criteria 


Xmin = Xj es. Xmax 
are placed in the appropriate class interval (cell) by the algorithm: 
Class Interval (Cell No.) = INT (12% — Xn) + 0.5| 


A data value falling outside the specified range is discarded when 
entered and this condition is indicated by a flashing-nines display. The 
maximum allowable data range is 10°? to 10°. 


Outputs of the program are: 
n = total number of x; entered within the range specified 


2X; 


m = mean value of the x, = = 


AZ 
o* = variance of the x, = t Ex! 7 Oe 


f; = frequency of class interval (cell) j = M 


where N; = number of items in cell j 
ee . Xmax — Xmin ); 
U; = upper limit of cell j = (Xnax Xan + Me, 
lf allentries occur in one cell, then the frequency is 1 and the cell number 
which is added to the frequency before display, will appear to be one 


greater than it should be; that is, if all observations fall in cell 7, the 
frequency of cell 7 will be displayed as 8.000. 


Example: Xmin = —6, Xmax = 6 
Vata values: —6,,--$, 4, —3, —2, —1,-0,1, 2,3. 4,5 
Number of entries= n= 12 (in the range Xin = X < Xmax) 
Mean = m =—0.5 
Variance = o7 = 11.91666667 
Frequency and upper limit versus cell number: 





Cell No. 1 2 3 4 2 6 
Frequency 

Upper Limit |—5.000 -—4.000 -—3.000 -—2.000 -—1.000 0.000 
Cell No. 4 8 9 10 11 12 


Frequency 





UpperLimit| 1.000 2.000 3.000 4.000 5.000 6.000 
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USER INSTRUCTIONS 


Enter lower limit 


Display variance 


Display frequency 


Display upper limit of 


cell of step 10 


Notes: 1) The format of the display in Step 10 is (cell number. frequency). The cell number is incremented by 1 with each press 
of CE] 
2) Step 11 must be performed immediately following step 10 for the cell in question. 
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PROGRAM LISTING 


29 *x! 


56 *rtn 
46 *LBL 46 *LBL 
14 D 

“st flg 

{ 

"LBL 









“Denotes 2nd function key 


REGISTERS 


a xa 


FLAGS 
SE tae EE CRE, RY, ee Se 
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PROGRAM LISTING 


86 *rset 
*LBL 





*Denotes 2nd function key 
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NORMAL DISTRIBUTION 
In order to calculate the standard normal distribution, the following 


probability function is used: 


{ 
V2Tr 

Then: Q(x) = 2(x)(bit + beat? + bat? + bat4 + bst5) + e(x) 
1 


e7X?*/2 





Z(x) = 


Where: t = 


1 + Ox 
le(x)| < 7.5 x 10-8 
p = .2316419 
bi = .319381530 
be = —.356563782 
b3 = 1.781477937 
b4 = —1.821255978 


bs = 1.330274429 


NOTES: 1. Z(—x) = Z(x), P(x) = 1 — Q(x) 
2,.\X|.<.21.0917 
3. Z(x) must be calculated before Q(x) 


Reference: Handbook of Mathematical Functions, National Bureau of 
Standards, 1964. 


Example: For x = 2.02, Z(x) = .0518635767, Q(x) = .0216916245 
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USER INSTRUCTIONS 
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PROGRAM LISTING 


46 *LBL 30 *Vx 65 43 
11 A 95 = 93 
22 INV 03 
50 *st fig 01 
O1 1 09 
80 *if pos 03 
ie 08 
94 +/- 01 
50 “st fig 
1 
"LBL 
ny) 








On wawod- We xX 

oo Coa oa © 0 Ga © A GO © 

h- ON WH? OWAaDND OO 
NWOONNFPHOAN? HH KWON NO 


© 
N 


“x o- & OO = @ h * 
ooemonocqcnan ods o 
OoO=- ABR ANHOWOAN DA ON 





tz 0 

113 
114 05 
113’ OS 
116 09 
The OF 
118 08 
119 85 
120 43 RCL 
121 00 
122 02 
123 45 
124 05 
125 65 
126 01 
127 93 
128 03 
129 03 
130 00 
131 02 
132 07 
138 04 
134 04 
135 02 
136 09 
137 95 = 

138 65 X 

139 43 RCL 


+ Omnogaan = 


™ RAO G 


~ 


On fF ANN OWWDe =| XK OA 


“Denotes 2nd function key 


142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
1359 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 


60 
01 
69 
96 
46 
69 
94 
85 
01 
95 
56 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


PROGRAM LISTING 


140 00 0 
141 01 1 
95 = 
“if flg 
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CHI-SQUARE DISTRIBUTION 


The chi-square density function is calculated by 


p,. f(x) = a 


where v = degrees of freedom, 0 < vp = 141,x > OandInx = gece 


The following series expansion calculates the cumulative distribution: 


os 
» 


P(x|v) = (5) (2) 1 i x resis ames] 
2 


Where: ['(v+1) = vI'(v) = v! = v(v—1)! 
and r(5) = ae 


NOTES: 1. If the degrees of freedom is entered as a non-integer, it is 

rounded up to an integer for the calculation. 

2. If v or x is negative or zero, a flashing 0 is returned. 

3. Av must be entered before an x, and both must be entered 
before P is calculated. 

4. A new x may be entered without changing v. 

9. An x must be entered each time P is to be calculated. 

V 


Example: Forv = 10, r(5) = 24 and if x = 8.5, f(x) = .0969533757, and 
P= .41988 16863 


Reference: An /ntroduction to Probability and Stochastic Processes; 
James L. Melsa, Andrew P. Sage; Prentice-Hall Electrical Eng. Series, 
1973 
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Enter Program (A and B) 
Enter degrees 
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PROGRAM LISTING 


000 46 *LBL 028 <5.4- 056 04 4 084 52 EE 


































































































































































001 43 RCL 029 48 *EXC |057 49 *PROD |085 57 “fix 
002 02 2 030 00 0 058 00 0 086 00 0 
003 44 SUM 031 06 6 059 03 3 087 52 EE 
004 00 0 032 95 = 060 46 *LBL 088 22 INV 
005 08 8 033 22 INV 061 42 STO 089 52 EE 
006 43 RCL 034 90 *ifzro |062 43 RCL 090 57 “fix 
007 00 0 035 43 RCL 063 00 0 091 09 9 
008 02 2 036 43 RCL 064 04 4 092 22 INV 
009 55 = 037 00 0 065 75 — 093 80 *ifpos 
010 43 RCL 038 06 6 066 93 094 10 *E’ 
011 00 0 039 49 *PROD |067 05 5 095 90 *ifzro 
012 08 8 040 00 0 068 95 = 096 10 *E’ 
013 65 X 041 05 5 069 22 INV 097 42 STO 
014 43 RCL 042 43 RCL 070 90 *ifzro |098 000 
015 00 0 043 00 0 071 44 SUM 099 07 7 
016 04 4 044 05 5 072 59 *z 100 55 = 
017 85 + 045 56 *rtn 073 30 *Vx 101 01 1 
018 42 STO 046 46 *LBL [074 49 *PROD |102 42 STO 
019 00 0 047 44 SUM 075 00 0 103 00 0 
020 04 4 048 43 RCL 076 03 3 104 03 3 
021 43 RCL 049 00 0 077 43 RCL 105 02 2 
022 00 0 050 04 4 078 00 0 106 95 = 
023 06 6 051 75 — 079 03 3 107 42 STO 
024 95 = 052 01 1 080 56 *rtn 108 00 0 
025 52 EE 053 95 = 081 46 *LBL 109 01 1 
026 22 INV 054 42 STO 082 42 


O2/ 52 EE 055 00 


*Denotes 2nd function key 
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PROGRAM LISTING 






















































112 04 4 140 56 * 

113) 73 - 141 46 a 
114 53 ( 142 12 B 
115 OF “Tix 143 22 INV 
116 00 O 144 52 EE 
tte (Oo = 145 22 INV 


118 93 
119 05 5 
120 54 ) 
121 52 EE 
122 22 INV 
123 52 EE 
124 57 *fix 
125 09 9 
126 95 = 
127 22 INV 
128 90 *ifzro 
129 42 STO 
130 43 RCL 
131 00 0 
132 01 1 
1o4 76 = 
134 01 1 
135 95 = 
136 29 *x! 
137 42 STO 
138 00 0 
139 


146 80 *ifpos 


“Denotes 2nd function key 
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t DISTRIBUTION 


The t distribution is evaluated by the following series expansion: 
: 0 + sin 6| cos O + 3 cos? ? +. + oe ae cos”- a] 


T 3 1-3... (v—2) 
for v>1 and odd 


A(t\v) = <P for v=1 


A(t]v) = 


A(t|v) = sin ol + Leos’ 64 ope! 64... 4 2 — cos’~? a | 


2 2:4 2°46 ... (v—2Z) 
for v even 
Where: v is the degrees of freedom 
t> 
t 
g= tan’ ' — = 
Vv 
Then: p= Alt +1 


NOTES: 1. Calculations must be done in radian mode. 
Example: For t= 1.341 and vy = 15, F = .9000625642 


Reference: Handbook of Mathematical Functions, Abramowitz and 
Stegun, National Bureau of Standards, 1964 
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ST1-17 





ST1-17 


PROGRAM LISTING 


02 2 


AY KEY 
46 *LBL 30 *Vx 
11 A 95 
42 92 
00 34 
01 42 
00 
03 
43 
00 
00 
55 - 
02 
75 
93 
05 
95 
20 
52 
57 
00 
95 
52 











PROGRAM LISTING 


ST1-17 





























113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
ras. 
126 
12/7 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 


“Denotes 2nd function key 


03 
she 
39 


112 00 O 


3 


* 
TT 


95 = 


10 
56 
46 
93 
01 
42 
00 
06 
o1 
42 
43 
00 
04 
10 
56 


*E’ 


140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
12 
153 
154 
155 
156 
be 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
































06 6 
55 + 
53 ( 
43 RCL 
00 0 
06 6 
85 + 
01 1 
) 
X 





























54 
65 
43 RCL 
00 0 

O3 3 

33 Cos 
40 *x? 
95 = 
42 STO 
00 0 

05 5 

44 SUM 
00 0 

04 4 

02 2 

44 SUM 
00 0 

06 6 

43 RCL 
00 0 
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F DISTRIBUTION 


The F distribution is most commonly used in terms of the ratio of two 
independent variables with chi-square distributions. 


If Y and Z are independent random variables distributed by chi-square 
Y/V ‘ss 
Z/v2 
called the variance ratio. Its distribution is the F distribution with 11 and 
v2 degrees of freedom. 





laws with v1 and v2 degrees of freedom respectively, then F = 


This program calculates Q(x), the probability of an event ‘x > X’’. 
Q(x) = Pr{x>xX} 
Q(x) is the “upper tail’ area on the F distribution curve. 


f(X,v1, 12) 





Q(0) = 1.0 (total area under curve = 1.0) 


This program calculates Q(x) using series expansions, provided that 
either 11 Or v2 is even. 


Conditions: x=0, 11 and v2 must be integers > 0 and one must be even. 
Q(x) is calculated according to the formulas: 





1. vi even 
Q(x) = Ke] ¥ 2 T 9.4 Fea, 
4 2 (vat+2) ... (vatrr—4)(1—K). * | 
2) 2 (Ie) 
2. v2 even 


vi(vit2) ... (v2tvi—4) K © | 
2:4... (v2—2) 


V1 2 
al vik |, vi(vit2) K jell 





-4- 


, 


Where: K = qo gates 
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ST1-18 


If both 11 and v2 are even, either formula may be used. Choosing the 
smaller of the two will reduce computation time, since fewer terms are 
generated by the series expansion. For example, if v1 = 4 and v2 = 16, 
use the ‘v1 even’’ key Ml to calculate Q(x) more quickly. 


Example: For v1 =5, v2=8: Q(1.0) = .4748134031, Q(9.5) = 0032423325 
For 11 = 8, v2 = 5:Q(1.0) = 5251865969 
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Si-ILS NOILNGIISIGS =am> | 
_| @A F DISTRIBUTION T1-18 
Ma FaE Se Tees: CS aE, aya 






STEP PROCEDURE 


NO | — 


co 


if vi iS even 
if v2 is even 


recall inputs 
of the inputs, re-enter 
the desired change 


(step 2) and continue 


to step 3. 


Compute Q(x): 
4 lf it is desired to 





ST1-18 


PROGRAM LISTING 


46 *LBL 42 

32 sin 
43 
00 
09 
85 
53 
24 
75 
43 
00 
06 
85 
02 
54 
90 







“Denotes 2nd function key 


02 x 
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PROGRAM LISTING 





“Denotes 2nd function key 
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BIVARIATE NORMAL DISTRIBUTION 


i gt a eter 
270.0, V1—p? s 2{1—p") 
Where: pu, = mean of x 
bly = mean of y 
ao, = Standard deviation of x (a, # 0) 
ao, = standard deviation of y (a, # 0) 
p =correlation coefficient between x and y (p? < 1) 
Error Indications: >, =0, 0, =0, p? = 1 


Reference: Handbook of Tables for Probability and Statistics, edited by 
W. H. Beyer, The Chemical Rubber Company, 1966 





a ait 


Example: pn, = 1,0, = 1.2 
by = 1.3, oy = .7 
p= 8 
f(2,1) = .0960208596 
f(1,1.3) = .2187813077 
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TEXAS INSTRUMENTS 


6L-LLS ISIC 1VINYON SLVINVAISG =a | 


USER INSTRUCTIONS 


[er Pran endy 
Enter o, 





DISPLAY 













Enter o,, 

8 
Enter x 

sz Enter y 

ok 


Calculate f(x,y) a ee 


— 
—_— 
< 
=< 
— 
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PROGRAM LISTING 


46 *LBL 056 
11 A 057 06 6 
STO 098 65 — 
059 43 
060 00 
061 00 
062 95 = 
063 55 + 
064 43 
065 00 
066 01 
067 95 
068 48 
069 00 
070 06 
071 56 
072 46 
073 17 
074 42 
075 00 
076 07 
Of? 73 
078 43 
079 00 
080 02 
081 95 = 
082 55 + 
083 43 RCL 







“Denotes 2nd function key 


REGISTERS 


13 
14 







01 oy 
ss ne I 7S Ee 
ol ES LS OM SR LE 
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PROGRAM LISTING 
































112 43 RCL 
113 00 O 
114 07 7 
115 40 *x? 
116 95 = 
117 55 + 
118 02 2 
119 55 + 
120 43 RCL 
121 00 0 
122 05 5 
123 95 = 
124 94 +/- 
125 22 INV 
126 23 Inx 
ler SO = 
128 43 RCL 
129 00 0 
130 05 5 
131 30 *Vx 
132 55 + 
las Gee 
134 55 + 
135 59 “9 
1386 55 + 
137 43 RCL 
138 00 0 
139 


“Denotes 2nd function key 
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LOGARITHMIC NORMAL DISTRIBUTION 


ee XV = ‘ey i miner. : 


Where: x is arandom variable whose logarithm is normally distributed 
with mean m and variance o?. (x > 0, a? # 0) 


For given m and o%, the following are computed: 


median = e™ 
mode = e‘™~?” 
mean, = ef? e 
variance = e +?™(e7— 1) 
f(x) 
Reference: Handbook of Probability and Statistics with Tables, 
Burlington and May, McGraw-Hill 

Example: m = 2, o?=1.5 


median = 7.389056099 f(.01) = .0000157407 
mode = 1.648721271 f(.5 ) = .058064287 
mean = 15.64263188 (1.0) = .0858628159 


variance = 851.9412262 
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CEE PE See ESS 







{| ererProvam tom | +t +i —s«dY 

2 | Computemedan fm —‘(cal] | 

3 | Computemade fo —‘(ce] | 

4 | Computemean | SS—«i CY 

"S| Compute varance «|_| vance 
LOAM <= WER ec 
Re i Ws: 
CACO SMS NE Bee 








STEP 


CS Se ermrater o 
| Renee or 
ite 






additional values of x 
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PROGRAM LISTING 


46 “LBL 95 = 
11 A 22 INV 
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WEIBULL DISTRIBUTION 


This program can be used to calculate: 
1. f(x1) = mnx™~ Ng(—mx) 
Where m>0, n>0, x>0 
2. Q(x1) =e™? 
3. X2 fora given Q, Q>0, such that 
Wiss (+ In oe 


m 
Example: m= .2,n=.9 
f(2) = .1156338734 
Q(2) = .6885184274 
if Q(x2) = .75, X2 = 1.497701748 












rig TEXAS INSTRUMENTS 


a a eee | 
|_| _#A®® WEIBULL DISTRIBUTION ST1-21 
Le ee a es Bee 


Q(x2) —> x2’ 


Enter m 


ESE 


4 [Compe ds 
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PROGRAM LISTING 


000 46 *LBL 084 07 7 














































































00 


001 11A 085 20 *1/x 
002 42 STO 086 95 = 
003 00 0 087 56 *rtn 
004 06 6 088 46 *LBL 
005 22 INV 089 87 *1’ 
006 80 *if pos 090 25 CLR 
O07 87 *1" 091 20 *1/x 
008 56 *rtn 092 56 *rtn 
009 46 *LBL 093 00 
010 12 B 094 00 
O11 42 STO 095 00 
012 00 0 096 00 
013 07 7 097 00 
014 22 INV 098 00 
015 80 “if pos 099 00 
016 87 *1’ 100 OO 
017 56 *rtn 101 00 
018 46 *LBL 102 00 
019 13 C 103 00 
020 42 STO 104 00 
021 00 0 105 00 
022 08 8 106 00 
023 51 SBR 107 00 
024 33 cos 108 OO 

65 109 00 

110 00 


S11-22 


POISSON DISTRIBUTION 


This program calculates the density function f(x) and the cumulative 
distribution P(x). 


m*e” ™ 
no 


Where:m>0O 
ie Oe 


Reference: Handbook of Tables for Probability and Statistics, edited by 
W. H. Beyer, The Chemical Rubber Company, 1966 


Example: m= 2.7 
f(0) = .0672055127 
f(8) = .0047075451 
P(8) = .9980863725 







TEXAS INSTRUMENTS 


we 


Se eee a 
| @A POISSON DISTRIBUTION ST1-22 
ree Se ea eS ae 





USER INSTRUCTIONS 


Enter Program (A only) 
Compute f(0) 


Compute P(x) 
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PROGRAM LISTING 


028 06 6 





029 
030 
031 
032 
033 
034 
035 
036 
037 
038 
039 
040 
041 
042 
043 
044 00 
045 03 
046 42 
047 00 


048 02 
049 00 
050 42 
051 00 
052 00 
053 46 
054 87 


AY 
00 
00 
73 
43 
00 
04 
35 
90 
32 
43 
00 
01 
99 
93 
43 
00 
00 
85 
01 
54 
65 
43 
00 
03 


— © 


95 = 


42 


©) 
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911-23 


BINOMIAL DISTRIBUTION 


This program evaluates the binomial density function defined by 


f(k;n,p) = (2) p* gq’ —” a oe 
O=<=p <1 
4 =41 +p) 
= 0 162.20. 


With mean: m = np and variance o? = npq 


Where: p = probability of success on each trial 
q = probability of failure 


For agiven k, n, and p, the function f(k;n,p) gives the probability of k 
successes and (n—k) failures, without regard to their order of 
occurrence, inn trials of an experiment. 
Results of the evaluation are: 

f(k;n,p) 

m 
ao? 


and the cumulative distribution 
k 
P(k)= 2 (%) p® gr-® 
k’=0 
Which is the probability of observing k or less successes inn trials, and 
P’(k) = = (Kk) p*qr-* 
k=k’ 


which gives the probability of observing at least k’ successes inn trials. 
For simplicity, f(k; n, p) is written as f(k). 
With k =0 
i{K)-= 10) = q? 
Fork >0 


Reference: An Introduction to Probability Theory and Its Applications, 
William Feller, John Wiley & Sons, 1968 


Example:n=7 p=0.45k=5 


m +=383.15 

o* = 41,7325 

(5) =.1172214914 
P(5) = .9642938078 
P’(5) = .1529276836 
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€Z-LLS NOILNGINLSIG IVINONIG =a | | 
#A« BINOMIAL DISTRIBUTION 
ee: Se ae 












ot | EierProoram hands 
aC es 
ra [ewes 
4 | compuiemean ie 
5 | Gomputevarance | 
oth 
icon 
ee 
neta 






é_[ Compute 

ef compact 
[Revers 7.8tor 
Pew vaeotk 


Note: All 9's (flashing) in the display following a data entry, indicates invalid data entry. 


| eleeleee 
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PROGRAM LISTING 


O1 1 
00 0 
46 *LBL 





ST1-23 


PROGRAM LISTING 


140 00 0 196 46 *LBL 
141'02 2 197 80 *ifpos 
142 94 198 00 O 
143 85 199 20 *1/x 
144 01 200 56 *rtn 
145 95 201 00 

146 42 202 00 

147 00 203 00 

148 03 204 00 

149 43 205 00 

150 00 206 00 

151 01 207 00 

152 65 208 00 

153 43 209 00 

154 00 210 

155 02 211 00 

156 95 = 212 00 

157 42 213 00 

158 00 214 00 

159 04 215 00 

160 56 216 00 

161 46 217 00 

162 14 = 218 00 

163 43 219 00 

164 00 220 00 

165 04 221 00 

166 65 222 00 

167 43 223 00 
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NEGATIVE BINOMIAL DISTRIBUTION 


This program evaluates the negative binomial density function defined by 


(ie) 

{ki fp) = k p'q* fe1i2...., 28 
Cag 1 
q =(1-— p) 
K=0,7,2.... 


With mean: m = and variance oa? = 8, 
Where: p = probability of success at each trial 
q = probability of failure 


For agiven k, r, and p, the function f(k; r,p) gives the probability that 
exactly k failures precede the rth success. 


Results of the evaluation are: 


f(k; r, p) 
m 
ao? 
and the cumulative distribution 
Pky = = (FET) pq 
k’=0 


which is the probability of observing k or less successes in r + k 
trials, and 


a oT oe Dh ae 
P= 2 (TT!) pa 


which gives the probability that there are at least k successes in 
r + K trials. 


For simplicity, f(k; r,p) is written as f(k). 
With k =0 


Fork > 0 
f(ky = VED Fh eK — 4) 
Reference: An /ntroduction to Probability Theory and Its Applications, 


William Feller, John Wiley & Sons, 1968 
Example:r=5 p=09 k=3 


m = .5555555556 P(3) = 0.99497565 
a2 = 6172839506 P’(3) = 0.0256915 
f(3) = 0.02066715 
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Panes SW Behe oe ae 


oC 





USER INSTRUCTIONS 


fer] vnoceoune [ewer [ness «omen 
Hi [tnerProoam areas [+t | ~<[ | 
LL ee RNS os as 
er Raa (2 A CN 














4 [conpueman fe ‘ey | [= 
"é | Compuetmy fk —ices| | ‘| 
compe | al ||| 
8 | comeueP [el | | 
9. 2 OU Lape ee ee 
ei aD SE AY TT 


Note: All 9’s (flashing) in the display following a data entry indicates invalid data entry. 
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PROGRAM LISTING 





911-24 


PROGRAM LISTING 


56 *rtn 
46 *LBL 
80 *if pos 
00 0 
20 "17x 
56 *rtn 
00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 





“Denotes 2nd function key 
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GEOMETRIC DISTRIBUTION 


For a geometric distribution with mean yp, this program will calculate: 
1. p, the probability of an event occurring 





2 the variance o? = __P 


oe 
3. f(x) = p(1 — p)* 
4. the cumulative distribution P(x) = 1— (1-—p)**! 
Where: 0 = p =1. 


Error Indications: A negative mean or x will cause the display to flash 
9.999999999 99 and those values will be ignored. 


Reference: Handbook of Probability and Statistics with Tables, 
Burington and May, McGraw-Hill 


Example:u=4.2 x=83 


p = .1923076923 
o* = 21.84 

(3) = .1013291726 

P(3) = .5744174749 
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Ceri ee ae dk) ©. 
_ [| #Am GEOMETRIC DISTRIBUTION ST1-25 





Enter Program (A only) 


Compute probability 


Compute f(x) 
Compute P(x P(x) 
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PROGRAM LISTING 
Oe de nT 


46 *LBL 

11 A 

22 INV 
“if pos 
i a 
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HYPERGEOMETRIC DISTRIBUTION 


This program evaluates the hypergeometric density function for 
given a, b, n: 


f(x) ae. 30) (nC —x) 
a+b@n 
Where: a, b, n are positive integers = 69 
a+b = 69 
X Sao=n +x =689 
x =O, 2... ..,62 
using the relation 
_ _ sip =1) 
(P+) oF hb—n+p +t 
p= 0, 1; 2/.5., 7) 


to find the cumulative distribution 
x—] 
POXy= 2 tp) 
p= 
The mean and variance are found using 
if 
a+b 


> Saba + Db — OB) 
~~ (a+ bila + 6 — 1) 





Example:a=10b=8n=7 


f(0) = P(0) = .0002513826 
f(4) = .3695324284 P(4) = .721719457 
f(7) = 0037707391 P(7) = 1. 

m = 3.888888889 

a? = 1118373275 
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9Z-11LS 1SIG DIYNLAWOSDYAdAH 2am sO 
_{._@Am™ HYPERGEOMETRIC DIST ST1-26 
Ceo ea ees Ba Be 


a,b n—> f(0) 









—s 


Enter a 
Enter b 


Compute mean 

For new x, goto 5. 

For new n, go to 4. 

For new a,b go to 2. 
0 | To compute m,o? for 


new a,b perform 


2,3, 7, and 8. 





ag 
az 
as 
9 
as 
i 
= 
Lad 


“Error indications after: 
Step 2 ifa<0, a>69, or non-integral 
Step 3 ifb<0, b>69, bnon-integral, or a+b>69 
Step 4 ifn<0, n>69, or non-integral 
step 5 if x<0, x>69, or non-integral 
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PROGRAM LISTING 


46 *LBL 03 3 20 + 084 16 *A’ 
Te *L’ 18 03 085 42 STO 
43 16 43 086 00 
00 43 00 087 O05 
00 00 02 088 01 
06 02 i = 089 94 
46 43 090 42 
18 00 091 00 
43 04 092 00 
00 04 093 43 
03 ga *x! 094 00 
96 95 095 08 
46 42 096 42 
11 097 00 
16 098 07 
42 099 46 
00 100 24 
01 101 42 
102 00 
103 06 
104 
105 
106 
107 
108 
109 
110 
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PROGRAM LISTING 


liz GS 5 196 01 1 
197 54 ) 
198 95 = 
199 56 *rtn 
200 46 *LBL 
201-17 "8 
202 43 RCL 
203 00 
204 01 1 
205 65 X 
206 43 RCL 
207 00 O 
208 04 4 
209 55 + 
210° 18 *C’ 
11 oS = 
2\2 56 *rtn 
213 00 0 
214 00 O 
215 00 O 
216 00 O 
217 00 0 
218 00 0 
0 
0 
0 
0 
0 












=) 






















219 00 
220 00 
221 00 
222 00 
223 








“Denotes 2nd function key 
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